CASE FILE 



ARR No. 3E28 



COPY 

NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



WARTIME REPORT 



ORIGINALLY ISSUED 

May 1943 as 
Advance Restricted Report 3E28 



THE EFFECT OF COMPRESSIBILITY ON TWO-DIMENSIONAL 
TUNNEL-WALL INTERFERENCE FOR A SYMMETRICAL AIRFOIL 

By Gerald E. Nitzberg 



Ames Aeronautical Laboratory 
Moffett Field, California 




WASHINGTON 



NACA WARTIME REPORTS are reprints of papers originally issued to provide rapid distribution of 
advance research results to an authorized group requiring them for the war effort. They were pre- 
viously held under a security status but are now unclassified. Some of these reports were not tech- 
nically edited. All have been reproduced without change in order to expedite general distribution. 



A-41 



NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS 



v ■ ADVANCE RESTRICTED REPORT 

TEE EFFECT OF COMPRESSIBILITY ON ? T WO-'DIMEN SI ONAL • 
TUNNEL- WALL INTERFERENCE FOR -A SYMMETRICAL AIRFOIL. 
■ . By Gerald E, ITit.zbe'rg/ 



SUMMARY 



The effective change in the velocity of flow past a 
wing section, caused by the presence of wind-tunnel walls, 
is known lor potential flow. This . theory is extended "by 
investigation. of the two-dimensional compressible flow 
past a thin JRankihe Oval. .It is shown that for a symmet- 
rical section at zero angle of attack the velocity incre- 
ment due to the tunnel walls in the incompressible case 
must be multiplied by the factor l/l-M 2 to. take account 
of compressibility effects. The Mach number, M, corre- 
sponds to conditions in the wind-tunnel test section with 
the model removed. 



' IKTBTODUCTIOIT 



Present day high-speed operation of wind tunnels 
renders it important to study the effects of compressibil- 
ity on the tunnel-wall interference. Several studies of 
the two-dimensional inviscid compressible flow in an in- 
finite stream over sections derivable from a circle have 
been made. Two recent papers have demonstrated a method 
for extending these studies to flow bounded by a channel. 
Ernest Lamia (reference l)- considered the incompressible 
potential flow past a circular cylinder in a channel as 
a first order approximation and used the Janz en-Ray ieigh 
method to find the compressible- inviscid flow past the 
circular cylinder in a channel. Von Hantzsche and Wendt 
(reference 2) have extended Lamia's study to the special 
family of ellipses in a channel having a maximum incom- 
pressible velocity past their surface ea/uai to twice the 
free stream velocity. 

These special cases are far removed from the usual 
conditions of wind-tunnel experimentation. A Rankine Oval 
having a thickness ratio comparable to a symmetrical 
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airfoil section gives a reasonable representation of the 
effects of profile thickness on the tunnel-wall velocity 
corrections for -an airfoil. A. Fage (reference 3) reports 
tne solution,by Sir Richard Giazebr ook, of the incompresfl- 
iole potential flow over a Rankine Oval in a channel. 
The present report takes this solution as a first order 
approximation and applies an analysis paralleling that of 
Lamia's to find the two-dimensional compressible . inyiscid 
flow past a thin Rankine Oval in a channel. 



THEORY 



In the Janzen-Rayle.igh method (reference l) ( the 
general stream function, and thus the velocity V at 
any point in the flow, is expressed as an infinite series 
function /of the- free stream Ma ah number M ( correspond- 
ing to the free stream velocity • V ),. as 



o 



so that 



where 




/9* A 2 /B\|; V 



The potential flow' stream function if satisfies 
Laplace's equation A \|/ Q = 0/ Introducing the conditions 
of continuous and: irrotat ional flow, leads .to" 

Therefore the second approximation to the general stream 
funct ion 

^ = ^ 0 + > 1 M*. , ' 
can be found from the. first approximation ii, = \!/ 0 . 



For adiabatic changes, to the order of M , 



and 



_P_ = 



7-1 ~ 



i . 



— 1 — K r 



7-1 = 1 - lM S 



consequent ly 



k A " 1 w 

.o . M 



?a.Si + 

P 2 

where 

p fluid density at any point in the stream 

p Q fluid density in the undisturbed stream 

p fluid density • at stagnation point 
s ■ 

The x velocity component at any point in the flow i 



•v = to a* j ^2 + [^o ^o_:_il + !*i 



P 3y 3y c dy 
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M 



when terms having powers of M greater than the squi 
are neglected. 



ANALYSIS 



The potential stream function for a Hankine Oval 
symmetrically placed in a channel (reference 3) is. 



- y + A tan 



si nil r sin m 



-cosh ?< - cosh r cos m J 



where 
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- t 



sinh r sin( *£) 



8 van- (_^=^-UlVE£_) 

: M - cosii r cos (~~y 

\ b / 



and all - angles must be taken to lie between o and tt. 
t thickness of the Rankihe" Oval 

s distance .between source and sink producing Rankine 
Oval 

V Q velocity of undisturbed stream' 

i = m = 5ny t r = HI 

■ H E E 

H breadth of channel ' 

x coordinate measured", along axis of the channel (major 
axis of oval) and having the origin midway between 
source and s ink ; . \ , J... - - — : - — • • ; ' 

y. coordinate measured perpendicular to the axis of the 
channel.- and having the origin midway- between the 
tunnel ytal is ■ 



Applying equation (l) gives 

— cosh I -~ 



a 2 i 2tt fE 3 " ■ 7 : 2E S , 

^ *i ~ "T - sxn m i ™ c °sk 1 ■( co .r .cosh I - cos m) 



where 



■-g- (,cosh t cosh r cos xn) sinn t + - cosh ^> 
P. 3) J 



•H 

••• 2 2 
D = (cosh 1 - cosh r cos m) + (sinh r sin m) 

To solve this partial differential equation by Lamia's 
method, it is necessary to obtain cos ..m as a finite 
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series in D, r, and I from the preceding equation. 
This imposes the limit (cosh 2 r - l) must be less than 1. 
For the sake of simplicity set cosii r.= i + e, where e 
is a small quantity. Consider a solution of equation (l) 
having the form 5 

*i = ™ h — -~— 7 

Then obtain the functions gp(V) 'by equating those terms 
in the two expressions for & Bn ty lm which have the same 
powers of D. If "terms in e 2 . are' neglected 

\|/. = I sin a 

2tt. ■ ■ 

where 

<1qE ^ d i cosh ■ I + (a/2) I sinh t (1 - cosh 2 I) 
I = „ r + " + 7571 

/ v 2 ed n E 2 \ ed.E 2 / cosh 2 l\ E 3 ^ 
v 2 2 / D V 3 ' 12 D 

■On the surface of the Eankine Oval >]> x = o, this provides 
the conditions necessary for evaluating the two constants 
do and d x . For convenience the coordinates for the ends 
of the Oval axes can be used 



= o when \ 



I = o, m = 



H 

■ l 



^ m = o, I = cosh X (-E + cosh r) 



This- completely determines and it is now possible to 

calculate the velocity at' any point in the flow from equa- 
tion (2) (and the analogous equation for the y velocity 
component ) . 

The characteristic point on the Rankine Oval may be 
taken to be the point x = o, y = t/2. For this point 

V ° x=o D :x=o x=o 

(3) 
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Lanla ? s result for a circle is the special case obtained 
from equation (3) by taking e = o. 



It can be shown that 



(4) 



(i": -'4} A-} ' ■ * (t/s) (5) 

The most frequent' applications have a thickness ratio, 
t/c, less- than 0.3, where. 



t 
c 



f/s. 



L . -1/ 



cot (t/s) 

therefore (4) would be negligible in comp.arison--w-tth- (~5~)~ 
. since a and- - b are of the same order. In such cases 
it is unnecessary to evaluate ^ x and . 

In reference 2 it is shown that' for an ellipse having t/c 
less than 0.3, in a channel with breadth so chosen as to 
give at x=o the incompressible velocity past the ellipse 

of twice the free stream velocity, ( — ) is very nearly 

s • 0 x=o 

2+3M . For 0,3 thickness ratio the velocity at x=o past 

a Rankine Oval is within 10 percent of that. past -an ellipse 
In equation (o) set M=o, and take ■ 7 ■ 

V o / z=o,M=o v ■ D2 J ... 



or 



D2 



then 



GfiL 



2+3M 



which gives a check on the validity of equation (6) , 

A further simplification in the application of equa- 
tion (S) is 



-2(t/s)[l+e/6(l+(t/ s ) 2 )] 



(7) 



(1+ (t/s) 8 ) \ _ -J-JLt ( i- 1( if c *o a A ] 



Comparison of this value with table I, reference 3, shows 
that for M=o» t/s less than 0.4, and s/E less than 0.5, 



the values for 



check within 0.1 percent. If an 



accuracy of 0.5 percent is sufficient, equation (7) can he 
reduced to 



B~2 



2_t 
"tTc 



1IC j 

~2t L 

TTC 



1 + 



TTC\ 2 

V Of/. 



12 



1 - 2t' 

TTC 



(s) 



Setting equation (8) in (6) it follows that at x^o 
the velocity past the Bankine Oval is 

.3 



c 



1 _ 2 1 



1 



12 ^1 - 2t./TTc //J 
/ TTC \ 



TIC 



+ m s j" — -iLns. — (i + Lki Ci + st^l \ 

Ll - 3t/ nc V 12 V TTC/^ J 



(9) 
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within the limits t/c less than 0.30 and c/h less than 
l/4. From equation (9) it is seen that the velocity past 



the Rankine Oval at x=o in an infinite stream 



(1 - •) 



is a function of Bach number, therefore in calculating the 
tunnel velocity correction it is necessary to consider the 
Mach number correction (reference l) 

(r r \ ., . 1°JL±L \ ♦ ( m + am) 2 (10) 

Writing equation (9) in the form of equation (10) 
gives as the tunnel velocity correction for thin bodies 

Atc \ , 

yav> = Cl + ^jSxi Y^-V (id 

Vv 0 y x=0 - 12 w/ 

In' reference 1, it is dpaonstrated that for a circle 

.3 



( 17 5. 

av\ , Cl + ii m 2 ] Li 
w x=0 * 12 



This/gives an indication oj| the magnitude of change in the 
multiplying constant of M for very thick bodies- . 

In reference 3 it is shown that for c/H less than 
0.5, the percentage tunnel velocity correction for all 
points on the body is nearly constant. If this limit is 
exceeded, the restriction of the tunnel walls causes a 
change in the effective shape of the body, because the 
velocity correction varies from point to point on the sur- 
face of the cylinder. 

A comparison of equation (11) with incompressible 
flow tunnel, velocity correction as given in reference 4 
shows that 

f*L\ ■ =(l + m s Y(4I) (12) 

^compressible ^ ' ^o incompressible 



The present development neglects powers of U greater 
than the square. 
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The mathematical difficulties of extending the pres- 
ent analysis to., higher Mach numbers are formidable. 
However 9 an investigation of the flow characteristics at 
Mach numbers approaching 1 suggests an extension of equa- 
tion (12). To this end, consider the two-dimensional 
flow past a thin symmetrical airfoil in free air. The 
conditions of conservation of mass.* and adiabatic com- 
pression lead to 



1 + 



M 



7-1 



1 + 

2 



7-1/ 



(13) 



where «' 0 and % are, respectively, the distances be- 
tween the same two stream lines ahead of the airfoil and 
at the position of maximum velocity over the airfoil. The 
velocities at these positions correspond to the Mach num- 
bers M 0 and M. 7 is the ratio of the specific heats. 

The maximum velocity past the surface of a thin sym- 
metrical section is of the same magnitude as the velocity 
ahead of the section. . Therefore, for small Mach numbers, 
equation ('13) indicates that the distances between stream- 
lines at these two positions are inversely proportional 
to their local Mach numbers. This corresponds .to a natural 
contraction of the streamlines. 

I 

As the Mach numbers approach 1, y- approaches 1» 
To demonstrate this» let o 



= 1 - 



and M = 1 - 6 , 



where e and 6 are small quantities, 
for air, equation (13) can be written 



Then, since 



7 = 1.4 



'o i - 8 L i + 



1 +0.2 (1 
0.2 (1 



_26) 
2€ N - J 



V x b J Vi - J 



The physical significance of this relation is that in the 
vicinity of the velocity of sound, the streamlines in 
passing from the undisturbed stream over the airfoil sec- 
tion would net contract in free air. The presence of 
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tunnel walls imposes a contraction, which alters the flow 
markedly. To the order of accuracy of this report, 1 + M 2 

is the series equivalent • t o 1 — - and the latter form 

1 - IT 

expresses the large magnitude of the. compressibility cor- 
rection when M approaches 1. Therefore, It seems advis- 
able to write equation (12) 

v j 2 Vv / 

0 compressible 1 - M 0 incompressible. 



This relation was derived for a thin Rankine Oval 
and indicates the order of correction for a thin symmet- 
rical airfoil section at zero angle of attack in a two- 
dimensional wind tunnel. 



Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moff ett .Field, Calif. ■ 
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